Flux-lattice melting in LaOi_a;Fa;FeAs: first-principles prediction 
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We report the theoretical study of the flux-lattice melting in the novel iron-based superconductor 
LaOo.tjFo.iFeAs and LaOo. 925^0. ovs-FeAs. Using the Hypernetted-Chain closure and an efliicient 
algorithm, we calculate the two-dimensional one-component plasma pair distribution functions, 
static structure factors and direct correlation functions at various temperatures. The Hansen- 
Verlet freezing criterion is shown to be valid for vortex-liquid freezing in type-II superconductors. 
Flux-lattice meting lines for LaOo.i)Fo.\FeAs and I/fflOo. 925^0. ovs^'eyls are predicted through the 
combination of the density functional theory and the mean-field substrate approach. 



PACS numbers: 74.25.Qt, 74.25.-q, 74.70.-b, 05.70.Fh 

I. INTRODUCTION 

It is well known that the practical application of su- 
perconductors crucially depends on the high electric- 
ity transport without dissipation, so the searching for 
new superconductors with a high transition tempera- 
ture has attracted considerable attention in material sci- 
ence. Recently, the novel superconductivity in iron-based 
layered superconductors has been reported, providing 
a new path to the high temperature superconductiv- 
ity. The LaOFeAs under doping with f ~ ions at the 
O^^ sites was first discovered to exhibit superconduc- 
tivity with Tc = 26K by Kamihara et al [H. Then the 
holes were introduced in LaOFeAs by partially substi- 
tuting La with S'r, the superconductivity was observed 
with Tc = 25 K 0. Shortly after these studies, 
was pushed up surprisingly to above 40 K when La in 
LaOi-xFxFeAs was substituted by other rare earth ele- 
ments 0, In later experiments, the maximum Tc at 
about 55 K was achieved in SmOi-xF^FeAs d, 0] and 
Gdi-xThxOFeAs Though these superconductors 
have a high critical temperature, the mechanism and the 
possible difference compared with copper-based high-Tc 
superconductors are still very controversial. Obviously, 
more experimental and theoretical work are needed to 
understand the type-II superconductivity as well as the 
microscopic mechanism in these new iron-based super- 
conductors. 

Since the vortex-lattice solid (glass) state without lin- 
ear resistivity is crucial for the application of high-Tc 
superconductors, the melting of the flux-lattice in bulk 
type-II superconductors is of fundamental and practical 

significance i, 13111111111,111, EE [13, Hill. The 

elastic theory combined with the phenomenological Lin- 
demann criterion has been used to study the flux-lattice 
melting in layered superconductors 20], which is based 
on the instability of the crystal. However, the regular 
flux-latticc structure loses simultaneously in the melt- 
ing, the information for vortex liquid is then lacking, so 
the Lindemann criterion has drawback in the study of 
the vortex melting. The density functional theory has 



been successfully used in the flux-lattice melting of lay- 
ered superconductors [9I [TtI. [TsI. IT9| . Due to the large 
anisotropy in these superconductors, weak Josephson in- 
teraction between layers can be neglected, both the in- 
terplane and intraplane interaction between the pancake 
vortices are electromagnetic in origin. The density func- 
tional method is usually well appropriate for layered su- 
perconductors with large anisotropy. Recently, a rather 
large anisotropy ~ 15 [2l| was found in LaOo,gFQ,iFeAs 
and LaOo. 925^0. 075^e^s (1111-type) by a bond structure 
calculation. The strong anisotropy in these systems was 
also confirmed recently in experiments . 

In this work, we study the flux-lattice melting in these 
two samples in the framework of the density functional 
theory and the mean-field substrate model proposed re- 
cently [23j . In our first-principles calculations, the input 
contains just a few parameters. The two-dimensional 
(2D) one-component plasma (OCP) pair distribution 
functions (PDF), 2D OCP static structure factors (SSF) 
and 2D OCP direct correlation functions (DCF) are cal- 
culated in the celebrated Hypernetted-Chain (HNC) clo- 
sure. The validity of the Hansen- Verlet freezing criterion 
is analyzed. The melting lines for LaOosFQ ^FeAs and 
LaOo. 925^0. 075^6^5 are plotted by using the mean-field 
substrate approach. 

The paper is organized as follows. Section II describes 
the calculation of 2D OCP PDF. In section HI, we ana- 
lyze the Hansen- Verlet freezing criterion and present the 
calculation of flux-lattice melting lines. A short summary 
is given in the last section. 



II. CALCULATION OF PAIR DISTRIBUTION 
FUNCTIONS 



The electromagnetic vortex-vortex interaction is given 
by a three-dimensional (3D) anisotropic pair potential. 
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In Fourier space, it reads [1, [T^ 



TXl,{T){kl + iA/d')sin\k,d/2) 
kl[l + klXUT) + {AXl,{T)/d^) sin2(fc,d/2] ' 



(1) 

where T = pd<j)l/47rXlf,{T) and /3 ^ 1/fcsT, kj_ and fc^, 
denote the component of k perpendicular and parallel 
to the c-axis respectively, Xab{T) = Xab{T — 0)/[l — 
(T/Tc)^]-^/^ is the planar penetration depth, (fy^ = hc/2e 
is the flux quantum and d is the layer spacing between 
FeAs layers. In our calculation, d is twice lattice param- 
eter c, i. e. d ~ 17. 5A[l|, [Hj. According to Eq. (1), 
the intralayer potential is repulsive while the interlayer 
potential is attractive. Some parameters measured in a 
muon spin relaxation study [24| are used in the present 
work , e. g. XabiT = 0) = 2540 A, = 26 K for 
LaOo.gFo.iFe^s and Xab{T = 0) = 3640 A, = 22 
K for LaOo. 925^0.075-^6^5. For convenience, the mean 
inter-particle spacing is taken as the length unit in the 
following. 

In the classical liquid theory, PDF is defined as [2^1 



g{r,r ) 



p{r, r) 



(2) 



where p(r) represents the probability to find a particle at 
r, p(r, r ) is the probability to find two particles at r and 
r respectively, and the pair correlation function h{r) is 
defined by h{r) = g{r) — 1. 

In the uniform and isotropic liquid, the Ornstein- 
Zernike relation can be written as [25j 

h{r) = Cir)+p J dr'C(|r-r'|)/i(r'), (3) 

while in the Fourier space, it becomes 

/i(k) = C(k)+pC(k)/i(k). (4) 

The HNC closure is widely used in the vortex system 
of type-II superconductors [3,[l3,[ll|- In this closure, the 
bridge function is set to zero and the DCF can be written 
as 



C(r) = exp[-l3V{r) + Y{r)] - Y{r) - 1, 



(5) 



where Y{r) = g{r) — C{r) — 1. Following the method 
proposed in Ref . [9| , we extract the correlation functions 
by splitting C(r) and Y{r) into a short- and a long-ranged 
part as 



c"(p) = cr(p) + Q"(p) 



r"(p) = n"(p) + i/"(p), 



(6) 



(7) 



here y,"(p) = (3V{p,nd), Cf(p) = -(3V{p,nd) and 
^sip) = C!sip)^n,o, P denotes the in-plane coordinate 
and n is the layer index In this way, the HNC equa- 
tion in the short range part can be written as 
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CO(p)=expK"(p)]-r,"(p)-l. 



(8) 
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FIG. 1: (a) 2D OOP PDF (b) 2D OCP SSF and (c) 2D 
OCP DCF for LaOi-^F^FeAs (x=0.1,0.075) at different tem- 
peratures. In all figures, from top to bottom at the first 
peaks: lOK for LaOo.gFo.iFeAs, UK for LaOo.gFo.iFeAs, 
lOK for LaOo.g25Fo.o75FeAs, 20K for LaOo.gFo.iFeAs, 14K 
for LaOo.92sFo.o75FeAs, and 20K for LaOo.925Fo.o75FeAs. 



The in-plane PDF is then given by 

g"{p) ^ CO(p) + r°(p)-l 

= cO(p) + nO(p)-i 

= exp{Y,'{p)). (9) 

Using an efficient algorithm proposed previously [11], 
we are able to calculate the short range part of 2D DCF 
and Y{p) successfully and then obtain the in-plane PDF. 
These functions are tabulated in a discrete form at about 
1200 points. Since the PDF decreases rapidly in the or- 
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FIG. 2; DCF based on the mean-field substrate model for 
LaOi-^R^FeAs (x=0.1,0.075) at T=14K, B=200Gs. 



der of the mean inter-particle spacing, they can be cut 
off in the real space to simplify the calculation. The cri- 
terion for this cutoff is that the converged functions are 
not influenced. In our calculation, all of the converged 
functions are stable and the relative difference is less than 
10~^. Actually, the value of PDF saturates at high field, 
consistent with the well-known fact that the vortex sys- 
tem decouples to a 2D system at high field. These PDF 
can be regarded as the 2D OCP one, which is shown in 
Fig. 1(a) for LaOi^^F^FeAs (x=0.1, 0.075). It displays 
an oscillating behavior and becomes very small at large 
P- 



III. FREEZING CRITERION AND MELTING 
LINES 

In the classical liquid theory, the SSF is defined as 

S{k) = l + p J dr{g{r) - l)exp{ik ■ r). (10) 

To provide a visual representation of the vortex struc- 
ture, we calculate the 2D OCP SSF according to Eq. 
(10). The 2D OCP DCF are then easily achieved with 
the relationship S{k) = 1/(1-C(fc)). Both 2D OCP SSF 
and 2D OCP DCF are shown in Figs. 1(b) and (c). In 
fact, the maximum values of these functions are found to 
decrease as temperature increases, shedding an insight of 
the 2D flux-lattice system in the melting process. Fur- 
thermore, all these functions for LaOo,QFo,iFeAs exhibit 
a larger oscillation than those for LaOo. 925^0.075-^6^5 at 
the same temperatures, indicating that the 2D lattice 
structure in LaO0.9F0.1Fe As is more clear than that in 
LaOo. 925 Fq. 075 Fe As. 

In the density functional approach first suggested by 
Ramakrishnan and Yussouff 27], the free energy function 
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FIG. 3; Flux-lattice melting lines for LaOi-xF^FeAs. Inset; 
Flux-lattice melting line for BSCCO calculated in Ref. ^18]. 

is given by 

l3/\n ^ [ dr[p{r)ln?^ - SpM] 
J Pi 

- IJdrJ dr'C{\r^r\) 

X [pir)~f^Mr)-f% (11) 

here p{r) is the density at r, J)^^ is the three-dimensional 
homogenous liquid density, and C(|r — r |) is the two- 
body fluid phase DCF. The two-body interactions be- 
tween the particles are expressed in the second integral 
on the right hand side. 

Recently, a mean-field substrate model was proposed 
to investigate the thermodynamic properties in layered 
superconductors [1^. Here we describe it briefly for 
completeness. Since the interlayer interaction between 
vortices is much weaker than the intralayer one in the 
strongly anisotropic superconductors, the vortices can be 
regarded as a collection of pancake vortices. In this way, 
the interplane and out-of-plane interaction could be sep- 
arated and the DCF can be written as [l3| 

C,ip)=dC'''ip)S{z)-^, (12) 

where Vz{p) is the out-of-plane interaction potential. In 
Fourier space, the effective DCF of the mean-field sub- 
strate model takes the following form 

Cr,(/0 = C^^{K) - J dz'P^. (13) 

We employ the above equation to extract the DCF 
based on the substrate model at different temperatures 
and magnetic fields. For instance, we present the DCF 
at T=14K and B=200Gs in Fig. 2. As proposed in 
Ref. [l3|, the single order parameter C^^(G) = 0.856 
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{{G = Sir'^p/Vsy^'^) can be used to determine the flux- 
lattice melting lines, which is equivalent to the nontriv- 
ial solution for fSAft ~ 0. On the other hand, the 
Hansen- Verlet freezing criterion is successfully used in 
the classical liquid theory [2^ . This criterion states that 
the liquid freezes when the first peak of the structure 
factor reaches a critical value. The density functional 
theory in the vortex system of type-II superconductors 
is also based on the classical liquid theory, thus it is 
necessary to analyze the validity of the Hansen- Verlet 
freezing criterion in this issue. Actually, the single or- 
der parameter C^^(G') = 0.856 is just relative to the 
2D Hansen- Verlet freezing criterion through the equation 
S{k) = 1/(1 — C(fc)). For the vortex system in layered su- 
perconductors with a large anisotropy, the liquid freezes 
when the structure factor at k = G reaches the value 
6.94. Interestingly, such a value is well consistent with 
the value of S(k) at low- field regime along the melting 
line obtained in previous work [18]. 

Finally, we employ the freezing criterion described 
above to locate the melting lines. Critical fields for 
LaOi-xFxFeAs {x = 0.1,0.075) at different tempera- 
tures are calculated, the melting lines are shown in Fig. 
(3). For comparison, the melting line for BSGGO ob- 
tained in Ref. [l^ is also given in the inset. One can 
observe that the monotonic meting lines for LaOFeAs 
are qualitatively alike but much lower than that for 
BSGGO. Recently, a quite large vortex hquid area 
in LaO{),<jFa.iFeAs compounds is observed in an ex- 
perimental study based on the measurement of Nernst 
signal [2^, consistent with the melting lines predicted 
in the present work. Interestingly, we find that in 
LaOFeAs samples the melting line is changed substan- 
tially with the variation of F~ content, the melting 
line for LaO{j,gFQ,iFeAs is much higher than that for 
LaOo,g2bFo.o75FeAs. 

The very recently discovered ternary iron-arsenides 
(Ba, K)Fe2S2 (122-type)t30| were suggested to be nearly 
isotropic [3l|. The present method to locate the melting 
lines of course can not be applicable to these supercon- 
ductors with low anisotropy. 

IV. SUMMARY 

The methodological difference between the present 
work and the previous study in Ref. [l3| lies in the deter- 



mination of 2D OCP direct correlation functions. In Ref. 
17] these functions were obtained from Monte-Carlo sim- 
ulations, while they are achieved ab inito from the micro- 
scopic vortex interaction in the present study. Therefore, 
the present approach is more intuitive, and the results are 
obtained without statistical errors like in Monte-Carlo 
simulations. The method adopted in the present work 
is also different from that in Ref. 9] where the density 
functional theory was first extended to the flux-lattice 
melting. A consistency to the Clausius-Clapeyron rela- 
tion [l7[ can be reached with the use of the mean-field 
substrate model in the present work, but can not be ob- 
tained in Ref. Q. 

We use the density functional theory and the mean- 
field substrate model to investigate the flux-lattice melt- 
ing in iron-based superconductors LaOi-xFxFeAs{x = 
0.1,0.075). The 2D OCP PDF, 2D OCP SSF and 2D 
OCP DCF for LaOo.gFo.iFeAs and LaOo.925^o.o75^'e^s 
are calculated at different temperatures via the HNC clo- 
sure and an efficient algorithm. It is shown that the melt- 
ing lines are quite low, well consistent with an experimen- 
tal study on the Nernst signals (29|. Furthermore, the 
melting line is observed to change considerably with the 
modulation of F^ content in i^-doped LaOFeAs com- 
pounds, indicating that the interaction between vortices 
crucially depends on the F^ content. More experimental 
studies are needed to test the present theoretical predi- 
cations. 



V. ACKOWNLEDGEMENTS 



This work was supported by National Natural Science 
Foundation of China under Grant Nos. 10574107 and 
10774128, PCSIRT (Grant No. IRT0754) in University 
in China, National Basic Research Program of China 
(Grant Nos. 2006CB601003 and 2009CB929104). 

t Corresponding author 



]1] Y. Kamihara, T. Watanabe, M. Hirano and H. Hosono, 

J. Am. Chem. Soc. 130, 3296 (2008). 
[2] H. H. Wen, G. Mu, L. Fang, H. Yang and X. Y. Zhu, 

Europhys. Lett. 82, 17009 (2008). 
[3] X. H. Chen, T. Wu, G. Wu, R. H. Liu, H. Chen and D. 

F. Fang, Nature 453,761 (2008). 
[4] G. F. Chen, Z. Li, D. Wu, G. Li, W. Z. Hu, J. Dong, 

P. Zheng, J. L. Luo, N. L. Wang, Phys. Rev. Lett. 100, 



247002 (2008). 

[5] Z. A. Ren, J. Yang, W. Lu, W. Yi, X. L. Shen, Z. C. Li, 

G. C. Che, X. L. Dong, L. L. Sun, F. Zhou and Z. X. 

Zhao, Europhys. Lett. 80, 57002 (2008). 
[6] Z. A. Ren, W. Lu, J. Yang, W. Yi, X. L. Shen, Z. C. 

Li, G. C. Che, X. L Dong, L. L. Sun, F. Zhou and Z. X. 

Zhao,Chin. Phys. Lett. 25, 2385 (2008). 
[7] R. H. Liu, G. Wu, T. Wu, D. F. Fang, H. Chen, S. Y. Li, 



5 



K. Liu, Y. L. Xie, X. F. Wang, R. L. Yang, C. He, D. L. 
Feng and X. H. Chen (2008), arXiv:0804.2105 
[8] C. Wang, L. J. Li, S. Chi, Z. W. Zhu, Z. Ren, Y. K. Li, 
Y. T. Wa X. Lin , Y. G. Luo, X. F. Xu, G. H. Cao and 
Z. A. Xu (2008). [aTXiv:0804.4290 i 
[9] S. Segupta, C. Dasgupta, H. R. Krishnamurthy, G. L 
Menon and T. V. Ramakrishman, Phys. Rev. Lett. 67, 
3444 (1991); G.I.Menon, C. Dasgupta, H. R. Krishna- 
murthy, T. V. Ramakrishman and S. Segupta, Phys. Rev. 
B 54, 16192 (1996). 

[10] E. M. Brandt , Phys. Rev. Lett. 63, 1106 (1989). 

[11] H. Safar, P. L. Gammel, D. A. Huse, D. J. Bishop, J. 
P. Rice and D. M. Ginsberg, Phys. Rev. Lett. 69, 824 
(1992). 

[12] W. K. Kwok, J. Fendrich, S. Fleshier, U. Welp, J. 
Downey and G. W. Crabtree, Phys. Rev. Lett. 72, 1092 
(1994).' 

[13] R. Cubitt, E. M. Forgan, G. Yang, S. L. Lee, D. Mck. 

Paul, H. A. Mool, M. yethiraj, P. H. Kes, T. W. Li, A. 

A. Menovsky, Z. Tarnawski and K. Mortensen, Nature 

365, 407 (1993). 
[14] C. Dasgupta and O. T. Vails, Phys. Rev. Lett. 91, 127002 

(2003); C. Dasgupta and O. T. Vails, Phys. Rev. B 74, 

184513 (2006); C. Dasgupta and O. T. Vails, Phys. Rev. 

B 76, 184509 (2007). 
[15] D. R. Nelson, Phys. Rev. Lett. 60, 1973 (1988). 
[16] A. Houghton, R. A. Pelcovits and A. Sudbo, Phys. Rev. 

B 40, 6763 (1989); X. Hu and Q. H. Chen, Phys. Rev. 

Lett. 92, 209701 (2004); Q. M. Nie, M. B. Luo, Q. H. 

Chen and X. Hu, Europhys. Lett. 71, 445 (2005). 
[17] A. De Col, G. I. Menon and G. Blatter, Phys. Rev. B 75, 

014518 (2007). 



[18] P. S. Cornagha and C. A. Balseiro, Phys. Rev. B 61, 784 
(2000). 

[19] X. Hu, M. B. Luo and Y. Q. Ma, Phys. Rev. B 72, 174503 
(2005). 

[20] A. Houghton, R. A. Pelcovits, and A. Sudbo, Phys. Rev. 

B 40, 6763 (1989). 
[21] D. J. Singh and M. H. Du, Phys. Rev. Lett. 100, 237003 

(2008). 

[22] G. F. Chen, Z. Li, G. Li, J. Zhou, D. Wu, J. Dong, W. 
Z. Hu, P. Zheng, Z. J. Chen, H. Q. Yuan, J. Singleton 
J. L. Luo and N. L. Wang, Phys. Rev. Lett. 101, 057007 
(2008). 

[23] M. J. W. Dodgson, A. E. Koshelev, V. B. Geshkenbein 
and G. Blatter, Phys. Rev. Lett. 84, 2698 (2000). 

[24] H. Luetkens, H. H. Klauss, R. Khasanov, A. Amato, R. 
Klingeler, I. Hellmann, N. Leps, A. Kondrat, C. Hess, 
A. Kohler, G. Behr, J. Werner and B. Biichner (2008), 
arXiv: 0804.3115. 

[25] J. P. Hansen and I. R. McDonald, Theory of Simple 
Liquids {3rd edition), (Academic Press, London, 2006). 

[26] M. J. Gillan, Molecular Physics 38, 1781 (1979). 

[27] T. V. Ramakrishman and M. Yussouff, Phys. Rev. B 19, 
2775 (1979). 

[28] J. P. Hansen and L. Verlet, Phys. Rev. 184, 151 (1969). 
[29] Z. W. Zhu, Z. A. Xu, X. Lin, G. H. Cao, C. M. Feng, G. 

F. Chen, Z. Li, J. L. Luo and N. L. Wang, New J. Phys. 

10, 063021 (2008). 
[30] M. Rotter, M. Tegel, D. Johrendt, (2008), arXiv: 

0805.4630. 

[31] H. Q. Yuan, J. Singleton, F. F. Balakirev, G. F. Chen, 
J. L. Luo, N. L. Wang (2008), arXiv: 0807.3137. 



